Unified description of correlations in double quantum dots 
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The two-level model for a double quantum dot coupled to two leads, which is ubiquitously used 
to describe charge oscillations, transmission-phase lapses and correlation-induced resonances, is 
considered in its general form. The model features arbitrary tunnelling matrix elements among the 
two levels and the leads and between the levels themselves (including the effect of Aharonov-Bohm 
fluxes), as well as inter-level repulsive interactions. We show that this model is exactly mapped onto 
a generalized Anderson model of a single impurity, where the electrons acquire a pseudo-spin degree 
of freedom, which is conserved by the tunnelling but not within the dot. Focusing on the local- 
moment regime where the dot is singly occupied, we show that the effective low-energy Hamiltonian 
is that of the anisotropic Kondo model in the presence of a tilted magnetic field. For moderate 
values of the (renormalized) field, the Bethe ansatz solution of the isotropic Kondo model allows us 
to derive accurate expressions for the dot occupation numbers, and henceforth its zero-temperature 
transmission. Our results are in excellent agreement with those obtained from the Bethe ansatz 
for the isotropic Anderson model, and with the functional and numerical renormalization-group 
calculations of Meden and Marquardt [Phys. Rev. Lett. 96, f4680f (2006)], which are valid for 
the general anisotropic case. In addition we present highly accurate estimates for the validity of the 
Schrieffer- Wolff transformation (which maps the Anderson Hamiltonian onto the low-energy Kondo 
model) at both the high- and low-magnetic field limits. Perhaps most importantly, we provide a 
single coherent picture for the host of phenomena to which this model has been applied. 
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I. INTRODUCTION 

The ongoing technological progress in the fabrication 
and control of nanoscale electronic circuits, such as quan- 
tum dots, has stimulated detailed studies of various 
quantum-impurity models, where a few local degrees of 
freedom are coupled to a continuum. Of particular in- 
terest are models with experimentally verifiable univer- 
sal properties. One of the best studied examples is the 
Anderson single impurity model^ which describes suc- 
cessfully electronic correlations in small quantum dots.^'^ 
The experimental control of most of the parameters of 
this model, e.g., the impurity energy level position or the 
level broadening due to hybridization with the contin- 
uum, allows for detailed investigations^ of the universal 
low-temperature behavior of the Anderson model. 

In this paper we study the low-energy behavior of a 
generic model, depicted in Fig. which pertains ei- 
ther to a single two-level quantum dot or to a double 
quantum dot where each dot harbors only a single level. 
The spin degeneracy of the electrons is assumed to be 
lifted by an external magnetic field. Several variants of 
this model have been studied intensely in recent years, in 
conjunction with a plethora of phenomena, such as many- 
body resonances in the spectral density,^ phase lapses 
in the transmission phase, ^'^ charge oscillations,^' and 
correlation-induced resonances in the conduct ance>ii*i^ 
Albeit being described by the same model, no clear link- 



age has been established between these seemingly differ- 
ent effects. The reason is in part due to the large num- 
ber of model parameters involved, which so far obscured 
a clear physical picture. While some exact statements 
can be made, these are restricted to certain solvable lim- 
itSfSi and are apparently nongeneric,^^ Here we construct 
a framework which encompasses all parameter regimes of 
the model, and enables a unified description of the var- 
ious phenomena alluded to above, exposing their com- 
mon physical origin. For the most interesting regime of 
strong fluctuations between the two levels, we are able to 
give: (i) explicit analytical conditions for the occurrence 
of transmission phase lapses; (ii) an explanation of the 
population inversion and the charge oscillationsSiiSiM (in- 
cluding a Kondo enhancement of the latter); (iii) a com- 
plete account of the correlation-induced resonanceaii as 
a disguised Kondo phenomenon. 

After introducing the details of the double-dot Hamil- 
tonian in Sec. Ill Al we begin our analysis by constructing 
a linear transformation of the dot operators, and a simul- 
taneous (generally different) linear transformation of the 
lead operators, such that the 2x2 tunnelling matrix be- 
tween the two levels on the dot and the leads becomes 
diagonal (with generally different eigenvalues). As a re- 
sult, the electrons acquire a pseudo-spin degree of free- 
dom which is conserved upon tunnelling between the dot 
and the continuum, as shown schematically in Fig. ^p. 
Concomitantly, the transformation generates a local Zee- 
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FIG. 1: A schematic representation of the double-dot system, along with its reduction in the local-moment regime to an effective 
Kondo model with a tilted magnetic field, (a) The model system: two localized levels coupled by tunnelling matrix elements 
to one another and to two separate leads. A constant magnetic flux induces phase factors on those elements. Spinless electrons 
residing on the two levels experience a repulsive interaction, (b) The mapping onto a spinful generalized Anderson model, with 
a tilted magnetic field and different tunnelling elements for spin-up and spin-down electrons, (c) The low-energy behavior of 
the generalized Anderson model is mapped onto an anisotropic Kondo model with a tilted magnetic field, htot ■ 



man magnetic field. In this way the original double-dot 
model system is transformed into a generalized Anderson 
impurity model in the presence of a (generally tilted) 
external magnetic field. This first stage is detailed in 
Sec. Ill Bl and AppendixEl 

We next analyze in Sec. Illll the low-energy properties 
of our generalized Anderson model. We confine ourselves 
to the local moment regime, in which there is a single 
electron on the impurity. The fluctuations of the pseudo- 
spin degree of freedom (which translate into charge fluc- 
tuations between the two localized levels in the original 
model) are determined by two competing effects: the po- 
larizing effect of the local magnetic field, and the Kondo 
screening by the itinerant electrons. In order to quanti- 
tatively analyze this competition, we derive an effective 
low-energy Kondo Hamiltonian, using Haldane's scaling 
procedurCf-^ together with the Schricffer- Wolff trans- 
formation and Anderson's poor man's scaling. This 
portion of the derivation resembles recent studies of the 
Kondo effect in the presence of ferromagnetic leads, al- 
though the physical context and implications are quite 
different. 

As is mentioned above, the tunnelling between the im- 
purity and the continuum in the generalized Anderson 
model is (pseudo) spin dependent. This asymmetry re- 
sults in two important effects: (a) different renormal- 
izations of the two local levels, which in turn generates 
an additional local magnetic fieldjii This field is not 
necessarily aligned with the original Zeeman field that 
is present in the generalized Anderson model, (b) An 
anisotropy of the exchange coupling between the con- 
duction electrons and the local moment in the Kondo 
Hamiltonian. However, since the scaling equations for 
the anisotropic Kondo modeliSiiS imply a flow towards 
the isotropic strong coupling fixed point, the low-energy 
behavior of the generalized Anderson model can be still 
described in terms of two competing energy scales, the 
Kondo temperature, T^, and the renormalized magnetic 
field, htot- Our two-stage mapping, double-dot ^ gener- 
alized Anderson model ^ anisotropic Kondo model (see 
Fig. ^ , allows us to obtain analytic expressions for the 
original model properties in terms of those of the Kondo 



model. We derive in Sec. llVl the occupation numbers on 
the two localized levels by employing the Bethe ansatz 
solution of the magnetization of a Kondo spin in a finite 
magnetic field>i2i2£ This solution also results in a highly 
accurate expression for the conductance based upon the 
Friedel-Langreth sum rule^^i Perhaps most importantly, 
it provides a single coherent picture for the host of phe- 
nomena to which our model has been applied. 

Examples of explicit results stemming from our general 
analysis are presented in Sec. First, we consider the 
case in which the tunnelling is isotropic, being the same 
for spin- up and spin-down electrons. Then the model is 
exactly solvable by direct application of the Bethe ansatz 
to the Anderson Hamiltonian. •^"■^^ We solve the result- 
ing equationsSSiS^ numerically and obtain the occupation 
numbers for arbitrary parameter values of the model, and 
in particular, for arbitrary values of the local Zeeman 
field. By comparing with the occupation numbers ob- 
tained in Sec. llVl from the Kondo version of the model, 
we are able to test the accuracy of the Schrieffer- Wolff 
mapping onto the Kondo Hamiltonian. We find that this 
mapping yields extremely precise results over the entire 
local-moment regime. This exactly solvable example has 
another virtue. It clearly demonstrates the competition 
between the Kondo screening of the local spin, which is 
governed by Tk, and the polarizing effect of the local field 
htot- This competition is reflected in the charging pro- 
cess of the quantum dot described by the original Hamil- 
tonian. We next proceed to apply our general method to 
the features for which the anisotropy in the tunnelling is 
relevant, notably the transmission phase lapses and the 
correlation-induced resonances^ In particular, we derive 
analytical expressions for the occupation numbers and 
the conductance employing the mapping onto the Kondo 
Hamiltonian. These analytical expressions give results 
which are in a very good agreement with the data pre- 
sented by Meden and Marquardt,^^ which was obtained 
by the functional and numerical renormalization-group 
methods applied to the original model. 

As our treatment makes extensive usage of the exact 
Bethe ansatz solutions for the impurity magnetization in 
the isotropic Kondo and Anderson models with a finite 
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magnetic field, all relevant details of the solutions are 
concisely gathered for convenience in Appendix FBI 



II. THE DOUBLE-DOT SYSTEM AS A 
GENERALIZED ANDERSON MODEL 

A. The model 

We consider spinless electrons in a system of two dis- 
tinct energy levels (a 'quantum dot'), labelled i = 1,2, 
which are connected by tunnelling to two leads, labelled 
a = L, R. This quantum dot is penetrated by a (con- 
stant) magnetic flux. The total Hamiltonian of the sys- 
tem reads 



Coulomb repulsion between electrons that occupy the 
two levels, eo ± A/2 are the (single-particle) energies on 
the levels, and 6/2 is the amplitude for tunnelling be- 
tween them. The phases tp^ and ipji, respectively, rep- 
resent the Aharonov-Bohm fluxes (measured in units of 
the flux quantum 2'nhc/e) in the left and in the right 
hopping loops, such that the total flux in the two loops 
\s LP = ipL + '^R [see Fig. [TJi] . 

Gauge invariance grants us the freedom to distribute 
the Aharonov-Bohm phases among the inter-dot coupling 
6 and the couplings between the dot levels and the leads. 
With the convention of Eq. 10}, the coupling between the 
quantum dot and the leads is described by the Hamilto- 
nian 



H. = 1-11 + Hd + Hid , 



(1) 



in which Tii is the Hamiltonian of the leads, Ti.d is the 
Hamiltonian of the isolated dot, and TCm describes the 
coupling between the dot and the leads. The system is 
portrayed schematically in Fig. 

Each of the leads is modelled by a continuum of nonin- 
teracting energy levels lying within a band of width 2D, 
with a constant density of states pm^ The corresponding 
Hamiltonian is given by 



kOL 



k^ka^ka ' 



(2) 



where cj,^ (Cka) creates (annihilates) an electron of wave 
vector k on lead a. The two leads are connected to two 
external reservoirs, held at the same temperature T and 
having different chemical potentials, fiL and iir, respec- 
tively. We take the limit f^L—^^J'R = in considering 
equilibrium properties and the linear conductance. 
The isolated dot is described by the Hamiltonian 



"Hd = [d\ dl] ■ £d ■ 



di 
d2 



+ U nin2 , 



where 
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Here, dj (dj creates (annihilates) an electron on the 
ith level, Ui = d\d^ are the occupation- number opera- 
tors (representing the local charge), U > denotes the 
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k 



di 
d2 



H.c. 



(5) 



where 



A = 



a-Ri aR2e 



aL2 

iip/2 



(6) 



Here the real (possibly negative) coefficients the 
tunnelling amplitudes for transferring an electron from 
the level i to lead a. Note that the Hamiltonian de- 
pends solely on the total Aharonov-Bohm flux ip when 
the interdot coupling b vanishes. Also, the tunnelling ma- 
trix A is assumed to be independent of the wave vector 
k. This assumption considerably simplifies the analysis 
while keeping the main physical picture intact. 



B. Mapping onto a generalized Anderson model 

The analysis of the model defined in Sec. IH Al emolovs 
an exact mapping of the Hamiltonian of Eq. onto a 
generalized Anderson Hamiltonian, which pertains to a 
single- level quantum dot, coupled to a spin-degenerate 
band of conduction electrons. We show in Appendix IXI 
that the model depicted in Fig. ^ is fully described by 
the Hamiltonian 
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Un^n^+Y,V„(^cl^d^ + R.. 



(7) 



schematically sketched Fig. which generalizes the lows for spin-dependent coupling between the dot and 
original Anderson modeli in two aspects. Firstly, it al- 
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the conduction band. A similar variant of the Anderson 
model has recently attracted much theoretical and ex- 
perimental attention in connection with the Kondo effect 
for ferromagnetic leads >i2iSS*2SiSL2& Secondly, it allows for 
a Zeeman field whose direction is inclined with respect 
to the "anisotropy" axis z. For spin-independent tun- 
nelling, one can easily realign the field along the z axis 
by a simple rotation of the different operators about the 
y axis. This is no longer the case once V\ 7^ Vj, which 
precludes the use of some of the exact results available 
for the Anderson model. As we show below, the main 
effect of spin-dependent tunnelling is to modify the effec- 
tive field seen by electrons on the dot, by renormalizing 
its z-component. 

The derivation of Eq. (Q is accomplished by a trans- 
formation known as the singular-value decomposition,^^ 
which allows one to express the tunnelling matrix A in 
the form 



A = R] 



Ft 

n 



(8) 



Here Ki and Rd are unitary 2x2 matrices, which are used 
to independently rotate the lead and the dot operators 
according to 



Rd 



dl 
d2 



Ri 



CkL 
CkR 



(9) 



To make contact with the conventional Anderson impu- 
rity model, we have labelled the linear combinations of 
the original operators [defined through Eqs. (|5J)] by the 
"spin" index a =] (+1) and a =[ (-1). 

The transformation (jSJ generalizes the one in which 
the same rotation R is applied to both the dot and the 
lead operators. It is needed in the present, more general, 
case since the matrix A generically lacks an orthogonal 
basis of eigenvectors. The matrices Rd and Ri can always 
be chosen uniquely (up to a common overall phase) such 
thalj^S (a) the tunnelling between the dot and the contin- 
uum is diagonal in the spin basis (so that the tunnelling 
conserves the spin); (b) the amplitudes V| > Vj > 
are real; and (c) the part of the Hamiltonian of Eq. {Tj) 
pertaining to the dot has only real matrix elements with 
hsmO > 0. The explicit expressions for the rotation ma- 
trices Rd and Ri as well as for the model parameters 
appearing in Eq. Q in terms of those of the original 
Hamiltonian are given in Appendix IXI 

It should be emphasized that partial transformations 
involving only one rotation matrix, either Rd or Ri , have 
previously been applied in this context (see, e.g., Refs. ,6 
and l3 j) . However, excluding special limits, both Rd and 
Ri are required to expose the formal connection to the 
Anderson model. A first step in this direction was re- 
cently taken by Golosov and GefenjS yet only on a re- 
stricted manifold for the tunnelling amplitudes aai- In 
the following section we discuss in detail the low-energy 
physics of the Hamiltonian of Eq. {Tj), focusing on the 
local-moment regime. Explicit results for the conduc- 



tance and the occupations of the levels are then presented 
in Sees. HVI and IVI 



III. THE LOCAL-MOMENT REGIME 

There are two limits where the model of Eq. has 
an exact solution:® (i) when the spin-down state is de- 
coupled in Eq. ITJ, i.e., when — hsinO — 0; (ii) when 
the coupling is isotropic, i.e., V| = V|. In the former 
case, ni is conserved. The Hilbert space separates then 
into two disconnected sectors with — and ni = 1. 
Within each sector, the Hamiltonian can be diagonalized 
independently as a single-particle problem. In the lat- 
ter case, one can always align the magnetic field h along 
the z axis by a simple rotation of the different operators 
about the y axis. The model of Eq. Q reduces then to 
a conventional Anderson model in a magnetic field, for 
which an exact Bethe ansatz solution is available, (This 
special case will be analyzed in great detail in Sec. lV Al ) 

In terms of the model parameters appearing in the orig- 
inal Hamiltonian, the condition V| = corresponds to 



\aLiaB.2\ = \aRiaL2\, and = /3 mod 2tt, 
whereas ~ Vi = V corresponds to 



(10) 



\a.Li\ ^ \aR2\, \aL2\^\aRi\, and = (tt + /?) mod 27r . 

(11) 

Here 



if aLiaL2aRiaR2 > 
TT if aLiaL2aRiaR2 < 



(12) 



records the combined signs of the four coefficients Uai^ 
Excluding the two cases mentioned above, no exact 
solutions to the Hamiltonian of Eq. ^ are known. Nev- 
ertheless, we shall argue below that the model displays 
generic low-energy physics in the "local-moment" regime, 
corresponding to the Kondo effect in a finite magnetic 
field. To this end we focus hereafter on T^,ri,h <C 
— eo, [/ + eo, and derive an effective low-energy Hamil- 
tonian for general couplings. Here T„ = npV^ is half the 
tunnelling rate between the spin state a and the leads. 



A. Effective low-energy Hamiltonian 

As is mentioned above, when V-\ = V| one is left with 
a conventional Kondo effect in the presence of a finite 
magnetic field. Asymmetry in the couplings, 7^ V|, 
changes this situation in three aspects. Firstly, the ef- 
fective magnetic field seen by electrons on the dot is 
modified, acquiring a renormalized 2r-component. Sec- 
ondly, the elimination of the charge fluctuations by 
means of a Schrieffer- Wolff transformation^^ results in 
an anisotropic spin-exchange interaction. Thirdly, a new 
interaction term is produced, coupling the spin and the 



5 



charge. Similar aspects have been previously discussed 
in the context of the Kondo effect in the presence of fer- 
romagnetic leads where the source of the asymmetry 
is the inequivalent density of states for conduction elec- 
trons with opposite spin.^^ Below we elaborate on the 
emergence of these features in the present case. 

Before turning to a detailed derivation of the effec- 
tive low-energy Hamiltonian, we briefly comment on the 
physical origin of the modified magnetic field. As is well 
known, the coupling to the continuum renormalizes the 
bare energy levels of the dot. For F-f , F^, h <C — cq, U + €o, 
these renormalizations can be accurately estimated using 
second-order perturbation theory in V^. For V-^ ^ V^, 
each of the bare levels €a ~ eo — \<^h cos 9 is shifted by a 
different amount, which acts in effect as an excess mag- 
netic field. Exphcitly, for T = and D » |eo|,C^ one 
obtainsiiii 



|eol 



(13) 



As eo is swept across — C//2, A/i^ cx F| — F; changes sign. 
Had |F| — F^ I exceeded h this would have dictated a sign- 
reversal of the z-component of the combined field as eo is 
tuned across the Coulomb-blockade valley. As originally 
noted by Silvestrov and Imry^i^ this simple but insightful 
observation underlies the population inversion discussed 
in Refs'o'To' ancT? for a singly occupied dot. We shall 
return to this important point in greater detail later on. 

A systematic derivation of the effective low-energy 
Hamiltonian for F|,F|,ft, ^ —eo,U + eo involves the 
combination of Anderson's poor-man's scaling^^ and the 
Schrieffer- Wolff transformationii^ For |eo| ^ U + eq, the 
elimination of high-energy excitations proceeds in three 
steps. First Haldane's perturbative scaling approacbi4 
is applied to progressively reduce the bandwidth from 
its bare value D down to Dsw kol ~ U + eo. Next 
a Schrieffer- Wolff transformation is carried out to elim- 
inate charge fluctuations on the dot. At the conclusion 
of this second step one is left with a generalized Kondo 
Hamiltonian [Eq. Ijltill below], featuring an anisotropic 
spin-exchange interaction and an additional interaction 
term that couples spin and charge. The Kondo Hamilto- 
nian also includes a finite magnetic field whose direction 
is inclined with respect to the anisotropy axis z. In the 
third and final stage, the Kondo Hamiltonian is treated 
using Anderson's poor-man's scalingi& to expose its low- 
energy physics. 

The above procedure is further complicated in the case 
where |eo| and U + eo are well separated in energy. This 
situation requires two distinct Schrieffer- Wolff transfor- 
mations: one at D^^ ~ niax{|eo|, U + eo} and the other 
at ~ min{|eo|,J7 + eo}. Reduction of the band- 

width from -Dg^ to i'g^" accomplished using yet 
another (third) segment of the perturbative scaling. It 
turns out that all possible orderings of |eo| and U + eo 
produce the same Kondo Hamiltonian, provided that F|, 
F| and h are sufficiently small. To keep the discussion as 
concise as possible, we therefore restrict the presentation 



to the case |eo| ^ U + eo- 

Consider first the energy window between D and Dsw^ 
which is treated using Haldane's perturbative scalingiii 
Suppose that the bandwidth has already been lowered 
from its initial value D to some value D' — De~^ with 
< I < ln{D/Dsw)- Further reducing the bandwidth 
to D'{1 — SI) produces a renormalization of each of the 
energies e|, e^, and U. Specifically, the z-component of 
the magnetic field, hz = ei — e^ , is found to obey the 
scaling equation 



dl 



Ft -F| 



1 



1-eJeo/D l + eJ{U + eo)/D 

(14) 

Here we have retained eo and U + eg in the denomina- 
tors, omitting corrections which are higher-order in F|, 
F|, and h (these include also the small renormalizations 
of Co- and U that are accumulated in the course of the scal- 
ing). The x-component of the field, hx — hsinO, remains 
unchanged throughout the procedure. Upon reaching 
D' — Dgw, the renormalized field hz becomes 



hi = hi 



Tt -Tx Dsw + U + eo 

— In — , (15) 

TT Dsw ~ eo 



where we have assumed D ^ |eo|: C^- 

Once the scale Dsw is reached, charge fluctuations on 
the dot are eliminated via a Schrieffer- Wolff transforma- 
tionjiS which generates among other terms also further 
renormalizations of e^. Neglecting h in the course of the 
transformation, one arrives at the following Kondo-type 
Hamiltonian, 

"Hk = ^ EfccL'^fccr + J-L {SxSx + SySy) + JzSzSz 



k.(7 



k,k',{T k,k',a 

hzSz — hxSx- 



k'a ■ 



(16) 



Here we have represented the local moment on the dot 
by the spin-^ operator 



(r being the Pauli matrices), while 



k.k' a, a' 



acr'^ka^k' a' 



(17) 



(18) 



are the local conduction-electron spin densities. The 
symbol ■■c\^Cf.,^: = cl^c^,^ - 6k^k'0{-ek) stands for nor- 
mal ordering with respect to the filled Fermi sea. The 
various couplings that appear in Eq. (|16|l are given by 
the explicit expressions 



pj± = 



1 1 

M U + eo 



(19) 
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pJz 



rt + Ti ( 1 



kol U + ea 



Ft -r. 



pVsc = 



1 



1 



C/ + eo 



rT±n 

47r 



1 1 

M " C^ + eo 



= h cos ( 



and 



hsinO. 



(20) 



(21) 



(22) 



(23) 



(24) 



Equations l|19|) - (|24|l are correct to leading order in r^, 
Fj^, and h, in accordance with the inequahty Fj, Fj, ft, <C 
|eo|,C/ + eo- In fact, additional terms are generated in 
Eq. (|16l) when h is kept in the course of the Schrieffer- 
Wolff transformation. However, the neglected terms 
are smaller than the ones retained by a factor of 
/i/ min{|eo|, J7 + eo} <C 1, and are not expected to al- 
ter the low-energy physics in any significant way. We 
also note that hz accurately reproduces the second-order 
correction to detailed in Eq. As emphasized 

above, the same effective Hamiltonian is obtained when 
|eo| and U + eg are well separated in energy, although 
the derivation is notably more cumbersome. In unify- 
ing the different possible orderings of |eo| and U + eg, 
the effective bandwidth in Eq. H16() must be taken to be 
Do ~ min{|eo|, ?7 + eo}- 



temperature. All relevant information on the bare spin- 
exchange anisotropy is contained for weak couplings in 
the microscopic form of Tk ■ 

The above picture is insensitive to the presence of 
weak potential scattering, which only slightly modifies 
the conduction-electron phase shift at the Fermi energy. 
As we show below, neither is it sensitive to the pres- 
ence of the weak couplings Ugc and u_ in Eq. IjlGI) . This 
observation is central to our discussion, as it enables a 
very accurate and complete description of the low-energy 
physics of Ti.K in terms of the conventional Kondo model 
in a finite magnetic field. Given the Kondo temperature 
Tk and the direction and magnitude of the renormalized 
field pertaining to Eq. H16|l . physical observables can be 
extracted from the exact Bethe ansatz solution of the 
conventional Kondo model. In this manner, one can ac- 
curately compute the conductance and the occupation of 
the levels, as demonstrated in Sees. II Vl andlvl 

To establish this important point, we apply poor- 
man's scalingi^ to the Hamiltonian of Eq. (|16|l . Of 
the different couplings that appear in Hk, only J^, J±, 
and hz are renormalized at second order. Converting 
to the dimensionless exchange couplings Jz = pJz and 
J± = pJ±-, these are found to obey the standard scaling 
equationsiSii^ 



dJz^ 
dl 
dJ± 
If 



JzJl. 



(25) 
(26) 



independent of v^c and v±. Indeed, the couplings ?;sc and 
v± do not affect the scaling trajectories in any way, other 
than through a small renormalization to hz '■ 



^ = ^0 e"' {^JzV- + 2iscV+^ 8 In 2. 



(27) 



B. Reduction to the Kondo effect in a finite 
magnetic field 

In addition to spin-exchange anisotropy and a tilted 
magnetic field, the Hamiltonian of Eq. (|16|l contains a 
new interaction term, Vsc, which couples spin and charge. 
It also includes spin-dependent potential scattering, rep- 
resented by the term v- above. As is well known, spin- 
exchange anisotropy is irrelevant for the conventional 
spin-i single-channel Kondo problem. As long as one 
lies within the confines of the antiferromagnctic domain, 
the system flows to the same strong-coupling fixed point 
no matter how large the exchange anisotropy is. SU(2) 
spin symmetry is thus restored at low energies. A fi- 
nite magnetic field h cuts off the flow to isotropic cou- 
plings, as does the temperature T. However, the resid- 
ual anisotropy is negligibly small if h, T and the bare 
couplings are small. That is, low-temperature thermo- 
dynamic and dynamic quantities follow a single generic 
dependence on T /Tk and h/T^, where Tk is the Kondo 



Here are the dimensionless couplings pw^ {p — sc, ±), 
and I equals \ti{Dq/ D') with D' the running bandwidth. 

As stated above, the scaling equations if^ - lf^ 
are identical to those obtained for the conventional 
anisotropic Kondo model. Hence, the Kondo couplings 
flow toward strong coupling along the same scaling tra- 
jectories and with the same Kondo temperature as in the 
absence of Vsc and v±. Straightforward integration of 
Eqs. yields 



Tk = Dq exp 



1 

'71 



■ tanh 



-1 e 



(28) 



Jj_. Here we have exploited the hierar- 



with C = 

chy > > in deriving Eq. (|28|l . In terms of the 
original model parameters appearing in Eq. {Tj) , Eq. H28|) 
takes the form 



Tk = Dq exp 



7reo(j7 + eo) I\ 
L2[/(Ft-Fx) ""f^ 



(29) 
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Equation (|29|l was obtained within second-order scaling, 
which is known to overestimate the pre-exponential fac- 
tor that enters Tk- We shaU not seek an improved ex- 
pression for Tff encompassing ah parameter regimes of 
Eq. 0. More accurate expressions will be given for the 
particular cases of interest, see Sec. fVl below. Much of 
our discussion will not depend, though, on the precise 
form of Tk- We shall only assume it to be sufhciently 
small such that the renormalized exchange couplings can 
be regarded isotropic starting at energies well above T]^. 

The other competing scale which enters the low-energy 
physics is the fully renormalized magnetic field: /itot = 



^tot ^ ~^ ^tot ^■ 



While the transverse field /if^j remains 



given by hsm6, the longitudinal field /i^^^ is obtained by 
integration of Eq. H27|l . subject to the initial condition 
of Eq. H23|) . Since the running coupling is a slowly 
varying function of I in the range where Eqs. 
apply, it can be replaced for all practical purposes by its 
bare value in Eq. H27|l . Straightforward integration of 
Eq. then yields 



Kot = hcos9 



U{U- 



2eo) 



(U + eo) 



2,2 ' 



(30) 



and 



where we have used Eqs. (|2Uf) - (|22|l for J^, Vsc, 
v±. Note that the third term on the right-hand side 
of Eq. is generally much smaller than the first two 
terms, and can typically be neglected. 

To conclude this section, we have shown that the 
Hamiltonian of Eq. {T)), and thus that of Eq. (QJ, is equiv- 
alent at sufficiently low temperature and fields to the 
ordinary isotropic Kondo model with a tilted magnetic 
field, provided that r|,r^ <^ |eo|,J7-|-eo. The relevant 
Kondo temperature is approximately given by Eq. (|29|l . 
while the components of htot = ^tot ^ + Mot ^ given 
by /iJqj = ft, sin 6* and Eq. ((201 . 



IV. PHYSICAL OBSERVABLES 

Having established the intimate connection between 
the generalized Anderson Hamiltonian, Eq. Q, and the 
standard Kondo model with a tilted magnetic field, we 
now employ well-known results of the latter model in or- 
der to obtain a unified picture for the conductance and 
the occupation of the levels of our original model, Eq. 
The analysis extends over a rather broad range of pa- 
rameters. For example, when U + 2eo — 0, then the 
sole requirement for the applicability of our results is for 
y/^ + W to be small. The tunnelling matrix A can be 
practically arbitrary as long as the system lies deep in the 
local-moment regime. The further one departs from the 
middle of the Coulomb-blockade valley the more restric- 
tive the condition on A becomes in order for fttot to stay 
small. Still, our approach is applicable over a surpris- 
ingly broad range of parameters, as demonstrated below. 



Unless stated otherwise, our discussion is restricted to 
zero temperature. 

A. Conductance 

At zero temperature, a local Fermi liquid is formed in 
the Kondo model. Only elastic scattering takes place at 
the Fermi energy, characterized by the scattering phase 
shifts for the two appropriate conduction-electron modes. 
For a finite magnetic field h in the z-direction, single- 
particle scattering is diagonal in the spin index. The 
corresponding phase shifts, S^{h) and Si{h), are given by 
the Friedel-Langreth sum rule^SLS Sa{h) = Ti{na), which 
when applied to the local-moment regime takes the form 



airMih) 



(31) 



Here M(h) is the spin magnetization, which reduceg^** in 
the scaling regime to a universal function of h/Tx, 



M{h)^MK{hlTK). 



(32) 



Thus, Eq. ^ becomes 5a{h) = 7r/2 + aTrMxih/TK), 
where Mxih/Tx) is given by Eq. 0\\ 

To apply these results to the problem at hand, one first 
needs to realign the tilted field along the z axis. This is 
achieved by a simple rotation of the different operators 
about the y axis. Writing the field /itot in the polar form 



htot = htot {sinOhX + cos 9hz) 
hsiiiO X + \ h cos 6 



TT |eo| 



(33) 

the lead and the dot operators are rotated according to 



CfeT 




CfcT 


= RhRi ■ 


CkL 











and 



h ■ 



— RhRd 



dl 

d2 



with 



cos{0h/2) sm{0h/2) 
-sin(0j2) cos(0j2) 



(34) 



(35) 



(36) 



Here Ri and Rd are the unitary matrices used in Eq. @ 
to independently rotate the lead and the dot operators. 
Note that since sin 6* > 0, the range of 6h is 9h G [0;7r]. 

The new dot and lead degrees of freedom have their 
spins aligned either parallel (d| and Cki) or antiparallel 
{dl and Cki) to the field /itot- In this basis the single- 
particle scattering matrix is diagonal. 



gi27rMK(/i,ot/TK) 







-i27rMjf (htot/Tjf) 



(37) 
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The conversion back to the original basis set of left- and 
right-lead electrons is straightforward, 



r t' 
t r' 



(38) 



providing us with the zero-temperature conductance G — 
(eV2^?i)|ip. 

Equations H37|l and H38|l were derived employing the 
mapping of Eq. |^ onto an effective isotropic Kondo 
model with a tilted magnetic field, in the Vsc,v± — > 
limit. Within this framework, Eqs. (|37|l and (|38|l are 
exact in the scaling regime, Tk/Dq ^ 1. The extent 
to which these equations are indeed valid can be appreci- 
ated by considering the special case hsinO = 0, for which 
there exists an exact (and independent) solution for the 
scattering matrix S in terms of the dot "magnetization" 
M = (n| - ni)/2 [see Eq. (@U below]. That solution, 
which is based on the Friedel-Langreth sum rulc^^ ap- 
plied directly to a spin-conserving Anderson model, re- 
produces Eqs. IpTzIl and in the Kondo regime. 



1. Zero Aharonov-Bohm fluxes 

Of particular interest is the case where no Aharonov- 
Bohm fluxes are present, where further analytic progress 
can be made. For LpL — ifiR = 0, the parame- 
ters that appear in the Hamiltonian of Eq. ^ are 
all real. Consequently, the rotation matrices Rd and 
Ri acquire the simplified forms given by Eqs. (|A29() 
and (|A32p (see Appendix ^ for details). Under 
these circumstances, the matrix product RuRi becomes 
±gir„(e^-HsHe,)/2gi,.T,(i-s„)/4^ j^^^ ^j^g elements of the 

scattering matrix [see Eq. (|38|l ] are 



t = t' = 



i sm[2TTMKihtot/TK)] sm{ei + SROh) , 

cos[2TTMK{htot/TK)] 

i sin[2TrMKihtot/TK)] cos{0i + sbA) ■ 

(39) 



Hence, the conductance is 



G = 



e2 

— sin'[2TTMK{htot/TK)] sin^Oi + srOh) , (40) 

where the sign sr and angle 0i are given by Eqs. (|X3l1) 
and (|A23p . respectively. All dependencies of the con- 
ductance on the original model parameters that enter 
Eq. Q are combined in Eq. (|40|) into two variables alone, 
01 + SR0h and the reduced field ft-tot/TV. In particular, 
6i is determined exclusively by the tunnelling matrix A, 
while sr depends additionally on the two dot parameters 
A and b. 

The conditions for a phase lapse to occur are particu- 
larly transparent from Eq. (|40ll . These lapses correspond 
to zeroes of t, and, in turn, of the conductance. There 
are two possibilities for G to vanish: either htot is zero, 
or 0i + SR9h equals an integer multiple of tt. For example. 



when the Hamiltonian of Eq. Q is invariant under the 
particle-hole transformation — > dj. and 



(which happens to be the case whenever V -|- 6^ = 
and U + 2eo = 0), then htot vanishes, and consequently 
the conductance vanishes as well. A detailed discussion 
of the ramifications of Eq. ()40|l is held in Sec. lVBl below. 



2. Parallel-field configuration 

For hsinO — 0, spin is conserved by the Hamiltonian 
of Eq. (TJ. We refer to this case as the "parallel-field" 
configuration, since the magnetic field is aligned with the 
anisotropy axis z. For a parallel field, one can easily gen- 
eralize the Friedel-Langreth sum rule^i to the Hamilto- 
nian of Eq. (|7|).^^ Apart from the need to consider each 
spin orientation separately, details of the derivation are 
identical to those for the ordinary Anderson model^i and 
so is the formal result for the T = scattering phase 
shift: 6cr — TrAiVcr, where ANg- is the number of dis- 
placed electrons in the spin channel a. In the wide-band 
limit, adopted throughout our discussion, AA^^- reduces 
to the occupancy of the corresponding dot level, (ria). 
The exact single-particle scattering matrix then becomes 



S 



J27tM 







,-i27rM 



Ri . 



(41) 



where M = {n^ — "■i)/2 is the dot "magnetization." 

Equation (|41|) is quite general. It covers all physical 
regimes of the dot, whether empty, singly occupied or 
doubly occupied, and extends to arbitrary fiuxes ipL and 
(pR. Although formally exact, it does not specify how 
the dot "magnetization" M and the total dot occupancy 
(7i| + n-i) relate to the microscopic model parameters 
that appear in Eq. ^ . Such information requires an ex- 
plicit solution for these quantities. In the Kondo regime 
considered above, (n| + rij) is reduced to one and M is 
replaced by ±Mi<-(/itot/Tif)- Here the sign depends on 
whether the field htot is parallel or antiparallel to the 
z axis (recall that htot > by definition). As a result, 
Eq. 1011) reproduces Eqs. (|S7|) - (|^ . 

To carry out the rotation in Eq. (|41|l . we rewrite it in 
the form 

S = e*''<"T+"i>i?t [cos(27rM) + i sin(27rM)T;,] Ri . (42) 

Using the general form of Eq. HA3|I for the rotation matrix 
Ri , the single-particle scattering matrix is written as S* = 
e*'^<"T+'H>5^ where 

S = cos(27rAf) -f-isin(27rM)cos6'/T2 

-I- i sin(27rM) sin 9i [cos (j)i + sin (j)i Ty] . (43) 

The zero-temperature conductance, G = (e^/27r?i)|i|^, 
takes then the exact form 



G = — sin2(27rM) sin^ 0i 



(44) 
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Two distinct properties of the conductance are appar- 
ent form Eq. H44I) . Firstly, G is bounded by sin^ 6i times 
the conductance quantum unit e^/27r?i. Unless 9i hap- 
pens to equal ±7r/2, the maximal conductance is smaller 
than /2'Kh. Secondly, G vanishes for M — and is 
maximal for M = ±1/4. Consequently, when M is tuned 
from Af « — 1/2 to M « 1/2 by varying an appropriate 
control parameter (for example, eg when 3> F|), then 
G is peaked at the points where M = ± 1/4. In the Kondo 
regime, when M — > ±M^(/itot/7V)j this condition is sat- 
isfied for /itot ~ 2.4Tx. As we show in Sec. lVB( this is the 
physical origin of the correlation-induced peaks reported 
by Meden and Marquardt.ii Note that for a given fixed 
tunnelling matrix A in the parallel-field configuration, 
the condition for a phase lapse to occur is simply for M 
to vanish. 



B. Occupation of the dot levels 

Similar to the zero-temperature conductance, one can 
exploit exact results of the standard Kondo model to ob- 
tain the occupation of the levels at low temperatures and 
fields. Defining the two reduced density matrices 



and 



Od 



{d{d,) (did,) 
{d\d,) (dld^) 



(dU^) (dU^) 

(d\d^) (dld^) 



these are related through 



Od - Rli^hOdRhRd ■ 



(45) 



(46) 



(47) 



Here RuRd is the overall rotation matrix pertaining to 
the dot degrees of freedom, see Eq. H35(l . 

At low temperatures, the mapping onto an isotropic 
Kondo model implies 



Od = 



where 



(nt> ' 
{ni) 



ntot/2 + 0-M, 



(48) 



(49) 



Here we have formally separated the occupancies {he) 
into the sum of a spin component and a charge compo- 
nent. The spin component involves the magnetization M 
along the direction of the total effective field /itot- The 
latter is well described by the universal magnetization 
curve Mxihtot/TK) of the Kondo model [see Eq. HBl|l ]. 
As for the total dot occupancy ntot, deep in the local- 
moment regime charge fluctuations are mostly quenched 
at low temperatures, resulting in the near integer valance 
ntot ~ 1. One can slightly improve on this estimate of 



ntot by resorting to first-order perturbation theory in Fo- 
(and zeroth order in h): 



«tot ~l + 



27r 



1 



1 



Co U + eo 



= 1 - 2pv+ . (50) 



This low-order process does not enter the Kondo effect, 
and is not contained in M^f (/itot/TV).^^ With the above 
approximations, the combination of Eqs. H47|l and H48|l 
yields a general formula for the reduced density matrix 



Od - ntot/2 + MK{htot/TK)RliRiT,RhR, 



Zero Aharonov-Bohm fluxes 



(51) 



As in the case of the conductance, Eq. (|51|) con- 
siderably simplifies in the absence of Aharonov-Bohm 
fluxes, when the combined rotation RuRd equals 
(sfi£0)^/2e'^«(^''+'"'^'^)''2e*'^^-(i-'"')/4 [see Eqs. ^ and 
29l]. Explicitly, Eq. if^ becomes 



7itot/2 + MK{htot/TK)cos{ed + seeh)T^ 

+ MK{htot/TK) sin(0<j + se0ft)T, , (52) 



Od = 



where the sign sg and angle 6d are given by Eqs. (jA30|) 
and HA23|I . respectively. 

Several observations are apparent from Eq. 152|1 . 
Firstly, when written in the original "spin" basis d\ 
and dl, the reduced density matrix Od contains the 
off-diagonal matrix element Mxihtot/TK) sm{6d + sgOh)- 
The latter reflects the fact that the original "spin" states 
are inclined with respect to the anisotropy axis dynam- 
ically selected by the system. Secondly, similar to the 
conductance of Eq. H40() . Od depends on two variables 
alone: Od + sgOh and the reduced field htot/Tx- Here, 
again, the angle 9d depends solely on the tunnelling ma- 
trix A, while the sign se depends additionally on A and b. 
Thirdly, the original levels d\ and dl have the occupation 
numbers 



m) = ntot/2 -I- MK{htot/TK) cos{0d ■ 
^ ntot/2 - MK{htot/TK) cos{0d 



(n2 



%)■ 



(53a) 
(53b) 



In particular, equal populations (ni) — (n2) are found 
if either /itot is zero or if Od -\- sgOd equals 7r/2 up to 
an integer multiple of tt. This provides one with a clear 
criterion for the occurrence of population inversion, ^494? ! 
i.e., the crossover from (ni) > (n2) to (n2) > (ni) or vice 
versa. 



2. Parallel-field configuration 

In the parallel-field configuration, the angle Oh is either 
zero or tt, depending on whether the magnetic field /itot is 
parallel or antiparallel to the z axis (recall that h sin = 
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htotsmOh = in this case). The occupancies (ni) and 
(7X2) acquire the exact representation 

{ni) = ntot/2 + M cos 9d , (54a) 
(ri2) = ritot/2 - M cos fld , (54b) 

where ntot is the exact total occupancy of the dot and 
M = (n| —ni)/2 is the dot "magnetization," defined and 
used previously (not to be confused with M = ±M). As 
with the conductance, Eqs. (|54|l encompass all regimes of 
the dot, and extend to arbitrary Aharonov-Bohm fluxes. 
They properly reduce to Eqs. (|53|l in the Kondo regime, 
when ntot « 1 [see Eq. (jHOU and M — > ±Mi^ (/^tot/7'i^)■ 
[Note that Eqs. (|53|l have been derived for zero Aharonov- 
Bohm fluxes.] 

One particularly revealing observation that follows 
from Eqs. H54(l concerns the connection between the phe- 
nomena of population inversion and phase lapses in the 
parallel-field configuration. For a given fixed tunnelling 
matrix A in the parallel-field configuration, the condi- 
tion for a population inversion to occur is identical to the 
condition for a phase lapse to occur. Both require that 
M = 0. Thus, these seemingly unrelated phenomena are 
synonymous in the parallel-field configuration. This is 
not generically the case when h^^^. ^ 0, as can be seen, 
for example, from Eqs. H40|l and (|53|) . In the absence of 
Aharonov-Bohm fluxes, the conductance is proportional 
to sin^ {61 + sjidh). It therefore vanishes for h^^^. ^ 
only if 9i + snOh = modTr. By contrast, the difference 
in populations (rii — 712) involves the unrelated factor 
cos{9d -\- sgOh), which generally does not vanish together 
with sin(0; + sj^Oh). 

Another useful result which applies to the parallel-field 
configuration is an exact expression for the T — con- 
ductance in terms of the population difference {rii — 71,2). 
It follows from Eqs. that M = {ni — / [2 cos 9 d) ■ 
Inserting this relation into Eq. 144() yields 

. 2 ( - ^^2) ^ . 2 a tKK\ 

G — — — sm sm 9i . (55) 

27rh \ cos9d J 

This expression will be used in Sec. for analyzing the 
conductance in the presence of isotropic couplings, and 
for the cases considered by Meden and Marquardtiii 



V. RESULTS 

Up until this point we have developed a general frame- 
work for describing the local-moment regime in terms 
of two competing energy scales, the Kondo temperature 
Tk and the renormalized magnetic field htot- We now 
turn to explicit calculations that exemplify these ideas. 
To this end, we begin in Sec. IV Al with the exactly solv- 
able case V| = Vi, which corresponds to the conventional 
Anderson model in a finite magnetic field. ^ Using the ex- 
act Bethe ansatz solution of the Anderson model,-" we 
present a detailed analysis of this special case with three 



objectives in mind: (i) to benchmark our general treat- 
ment against rigorous results; (ii) to follow in great detail 
the delicate interplay between the two competing energy 
scales that govern the low-energy physics; (iii) to set the 
stage for the complete explanation of the charge oscilla- 
^jpjjg9ji0ji3 g^j-jj ^]-^g correlation-induced resonances in the 
conductance of this deviceiii^i^ 

We then proceed in Sec. IV Bi to the generic anisotropic 
case V| ^ V|. Here a coherent explanation is pro- 
vided for the ubiquitous phase lapses,® population in- 
versioufSiiS, and correlation-induced resonance o^^i^^ that 
were reported recently in various studies of two-level 
quantum dots. In particular, we expose the latter res- 
onances as a disguised Kondo phenomenon. The general 
formulae of Sec. IIVI are quantitatively compared to the 
numerical results of Refflll The detailed agreement that 
is obtained nicely illustrates the power of the analytical 
approach put forward in this paper. 



A. Exact treatment of Vi- = V| 

As emphasized in Sec. IIIII all tunnelling matrices A 
which satisfy Eq. (|11|) give rise to equal amplitudes 
= Vi = V within the Anderson Hamiltonian descrip- 
tion of Eq. Q . Given this extra symmetry, one can al- 
ways choose the unitary matrices Ri and Rd in such a 
way that the magnetic field h points along the z direc- 
tion [namely, cos9 = 1 in Eq. Q]. Perhaps the simplest 
member in this class of tunnelling matrices is the case 
where uli = -aL2 = clri = ai?2 = V/V2, tpL = = 
and 6 = 0. One can simply convert the conduction- 
electron operators to even and odd combinations of the 
two leads, corresponding to choosing 9i — Tr/2 + 9d- De- 
pending on the sign of A, the angle 9d is either zero (for 
A < 0) or TT (for A > 0), which leaves us with a conven- 
tional Anderson impurity in the presence of the magnetic 
field h — |A| z. All other rotation angle that appear in 
Eqs. (|A2|) and ljA2p (i.e., x's and 0's) are equal to zero. 
For concreteness we shall focus hereafter on this particu- 
lar case, which represents, up to a simple rotation of the 
dl and c^^ operators, all tunnelling matrices A in this 
category of interest. Our discussion is restricted to zero 
temperature. 



1. Impurity magnetization 

We have solved the exact Bethe anstaz equations nu- 
merically using the procedure outlined in Appendix ^ 
Our results for the occupation numbers (n^) and the 
magnetization M — {nij—n^)/2 are summarized in 
Figs. 121 and 13 Figure [21 shows the magnetization of the 
Anderson impurity as a function of the (average) level 
position eg in a constant magnetic field, h = A = 10~^U. 
The complementary regime eo < —U/2 is obtained by 
a simple reffection about cq = ~U/2, as follows from 
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FIG. 2: (Color online) Magnetization of the isotropic case as a 
function of eo: exact Bethe ansatz curve and comparison with 
different approximation schemes. Black symbols show the 
magnetization M derived from the exact Bethe ansatz equa- 
tions; the dashed (red) line marks the result of first-order per- 
turbation theory in F (Ref0, divergent at eo = 0); the thick 
(blue) line is the analytical formula for the magnetization in 
the Kondo limit, Eq. jnU, with Tk given by Eq. The 
model parameters are F/f/ = 0.05, A/U = 10"^ and T = 0. 
The upper right inset shows the same data but on a linear 
scale. The lower left inset shows the magnetization M as a 
function of the magnetic field h — A at fixed eo/U — —0.2. 
The universal magnetization curve of the Kondo model well 
describes the exact magnetization up to M « 0.42 (lower 
fields not shown), while first-order perturbation theory in F 
fails from M ~ 0.46 downwards. 



the particle-hole transformation do- —f and — > 

— cj,_^. The Bethe ansatz curve accurately crosses over 
from the perturbative domain at large eo S> F (when 
the dot is almost empty) to the local-moment regime 
with a fully pronounced Kondo effect (when the dot is 
singly occupied). In the latter regime, we find excel- 
lent agreement with the analytical magnetization curve 
of the Kondo model, Eq. (|B1|) , both as a function of eo 
and as a function of the magnetic field A (lower left in- 
set to Fig. 121) . The agreement with the universal Kondo 
curve is in fact quite surprising in that it extends nearly 
into the mixed-valent regime. As a function of field, the 
Kondo curve of Eq. (|Bip applies up to fields of the order 
of VTU ^Tk. 



2. Occupation numbers and charge oscillations 

Figure 13 displays the individual occupation numbers 
(ni) and (n2) as a function of eo, for a series of con- 
stant fields h — A. In going from large eo ^ F to large 
— (eo + U) ^ F, the total charge of the quantum dot 
increases monotonically from nearly zero to nearly two. 
However, the partial occupancies (ni) and (712) display 
nonmonotonicities, which have drawn considerable theo- 
retical attention lately AiSii^ As seen in Fig. O the non- 
monotonicities can be quite large, although no popula- 



FIG. 3: (Color online) The occupation numbers (ni) [solid 
(blue) lines] and (712) [dotted (red) lines] versus eo, as obtained 
from the solution of the exact Bethe ansatz equations. In go- 
ing from the inner-most to the outer-most pairs of curves, the 
magnetic field h = A increases by a factor of 10 between each 
successive pair of curves, with the inner-most (outer-most) 
curves corresponding to A/U = 10~^ i^/U = 0.1). The re- 
maining model parameters are F/U — 0.05 and T = 0. Non- 
monotonicities are seen in the process of charging. These are 
most pronounced for intermediate values of the field. The evo- 
lution of the nonmonotonicities with increasing field is tracked 
by arrows. The dashed black lines show the approximate val- 
ues calculated from Eqs. 1531 and I50II based on the mapping 
onto the Kondo Hamiltonian (here 9h = and 9d = tt). 



tion inversion occurs for F^ = T^. 

Our general discussion in Sec. lIIII makes it is easy to in- 
terpret these features of the partial occupancies (rii). In- 
deed, as illustrated in Fig.O there is excellent agreement 
in the local-moment regime between the exact Bethe 
ansatz results and the curves obtained from Eqs. I|53() 
and H5U|I based on the mapping onto the Kondo Hamil- 
tonian. We therefore utilize Eqs. (|53|l for analyzing the 
data. To begin with we note that, for F| = F|, there is no 
renormalization of the effective magnetic field. The lat- 
ter remains constant and equal to ft. = A independent of 
cq. Combined with the fact that cos{9d + sg9h) = — 1 
in Eqs. H53I) . the magnetization M — (ri|— n|)/2 = 
(n2 — ni)/2 depends exclusively on the ratio A/Tk- The 
sole dependence on eo enters through Tk, which varies 
according to Eq. (jB4|) . Thus, M is positive for all gate 
voltages eo, excluding the possibility of a population in- 
version. 

The nonmonotonicities in the individual occupancies 
stem from the explicit dependence of Tk on the gate 
voltage eo. According to Eq. (|B4p . Tk is minimal in the 
middle of the Coulomb-blockade valley, increasing mono- 
tonically as a function of |eo + U/2\. Thus, A/Tk, and 
consequently M, is maximal for eo ~ — C//2, decreasing 
monotonically the farther eo departs from —U/2. Since 
ntot ~ 1 is nearly a constant in the local-moment regime, 
this implies the following evolution of the partial occu- 
pancies: (rii) decreases ((712) increases) as eo is lowered 
from roughly zero to —U/2. It then increases (decreases) 
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FIG. 4: (Color online) The exact conductance G [in units 
of eV(27rft)] versus eo, as obtained from the Bethe ansatz 
magnetization M and Eq. 1551 with 9i = 3n/2 and 6d ~ it. 
Here A/U equals 10"^ [full (black) line], lO"* [dotted (red) 
line], 10""^ [dashed (green) line] and 0.1 [dot-dashed (blue) 
line]. The remaining model parameters are T/U = 0.05 and 
T = 0. Once A exceeds the critical field « 2.4rip'", the 
single peak at eo = — is split into two correlation-induced 
peaks, which cross over to Coulomb-blockade peaks at large 
A. 



as eo is further lowered toward ^U. Combined with the 
crossovers to the empty-impurity and doubly occupied 
regimes, this generates a local maximum (minimum) in 
(m) (('^2)) near eo (eo -U). 

Note that the local extremum in (n^) is most pro- 
nounced for intermediate values of the field A. This can 
be understood by examining the two most relevant en- 
ergy scales in the problem, namely, the minimal Kondo 



temperature T^'" = Tj 



-U/ 



2 and the hybridization 



width r. These two energies govern the spin susceptibil- 
ity of the impurity in the middle of the Coulomb-blockade 
valley (when sq = —U /2) and in the mixed- valent regime 
(when either eo ~ or e « — ?7), respectively. The charg- 
ing curves of Fig. |31stem from an interplay of the three 
energy scales A, Tjp'" and F as described below. 

When A <C Tj^'", exemplified by the pair of curves cor- 
responding to the smallest field A = lO'^U « 0.24T]^^"' 
in Fig. |31 the magnetic field remains small throughout 
the Coulomb-blockade valley and no significant magne- 
tization develops. The two levels are roughly equally 
populated, showing a plateaux at (ni) « (712) ~ 1/2 
in the regime where the dot is singly occupied. As A 
grows and approaches TjJ?™, the field becomes sufficiently 
strong to significantly polarize the impurity in the vicin- 
ity of eo = — J7/2. A gap then rapidly develops between 
(ni) and (712) near eo — —U/2 as A is increased. Once 
A reaches the regime TjJ?'" <C A ^ F, a crossover from 
h ^ Tk (fully polarized impurity) to h <^ Tk (unpo- 
larized impurity) occurs as eg is tuned away from the 
middle of the Coulomb-blockade valley. This leads to 
the development of a pronounced maximum (minimum) 
in (ni) ((n2)), as marked by the arrows in Fig. 13 Fi- 
nally, when /i > F, the field is sufficiently large to keep 



FIG. 5: (Color online) The exact occupation numbers {rn) 
and conductance G [in units of e^/(27rft)] as a function of £2, 
for T = 0, r/17 = 0.2 and fixed ei/f/ = -1/2. The population 
inversion at £2 = ei leads to a sharp transmission zero (phase 
lapse). Note the general resemblance between the functional 
dependence of G on €2 and the correlation-induced resonances 
reported by Meden and Marquardliii for Fj / Fx (see Fig.|SJ. 



the dot polarized throughout the local-moment regime. 
The extremum in (n;) degenerates into a small bump in 
the vicinity of either eo ~ or eg ~ — C/, which is the non- 
monotonic feature first discussed in Ref^. This regime 
is exemplified by the pair of curves corresponding to the 
largest field A = O.lf/ = 2F in Fig. O whose parameters 
match those used in Fig. 2 of Rcf 9. Note, however, that 
the perturbative calculations of Ref 9 will inevitably miss 
the regime T]J?'" ^ A ^ F where this feature is large 



3. Conductance 

The data of Fig. O can easily be converted to conduc- 
tance curves by using the exact formula of Eq. I|55|) with 
Oi — 37r/2 and Od = tt. The outcome is presented in 
Fig. El The evolution of G(eo) with increasing A is quite 
dramatic. When A is small, the conductance is likewise 
small with a shallow peak at eo = —U 12. This peak 
steadily grows with increasing A until reaching the uni- 
tary limit, at which point it is split in two. Upon further 
increasing A, the two split peaks gradually depart, ap- 
proaching the peak positions eo « and eo ~ — C/ for 
large A. The conductance at each of the two maxima 
remains pinned at all stages at the unitary limit. 

These features of the conductance can be naturally un- 
derstood based on Eqs. ^ and lO- When A < T]^"\ 
the magnetization M « A/(27rTff) and the conductance 



G ~ (A/Tff) e /(27r?i) are uniformly small, with a peak 
at eo = —U/2 where Tk is the smallest. The conductance 
monotonically grows with increasing A until reaching the 



critical field A 



2.4TJI?'", 



where M|f„^_ 



■u/2 



1/4 



and G'|co=-t//2 = e?/{2'n:h). Upon further increasing A, 
the magnetization at eo = — C//2 exceeds 1/4, and the as- 
sociated conductance decreases. The unitarity condition 
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M = 1/4 is satisfied at two gate voltages ej^^x symmetric 



about — J7/2, defined by the relation 
Eq. l)B4p one obtains 



A/2.4. From 




r2 



2TU 



In 



ttA \ 



2.4V2rc7, 



(56) 



The width of the two conductance peaks, Ae, can be 
<C A <C r from the inverse of the 



estimated for Tf"" 



derivative d{A/TK)/deo, evaluated at eg 



Ae 



TU 



TrlCmax 



U/2\ 



It yields 



(57) 



Finally, when A > F, the magnetization exceeds 1/4 
throughout the local-moment regime. The resonance 
condition M = 1/4 is met only as charge fluctuations be- 
come strong, namely, for either ~ or eo ~ — ?7. The 
resonance width Ae evolves continuously in this limit to 
the standard result for the Coulomb-blockade resonances, 
Ae-F. 

Up until now the energy difference A was kept constant 
while tuning the average level position cq. This proto- 
col, which precludes population inversion as a function 
of the control parameter, best suits a single-dot realiza- 
tion of our model, where both levels can be uniformly 
tuned using a single gate voltage. In the alternative re- 
alization of two spatially separated quantum dots, each 
controlled by its own separate gate voltage, one could 
fix the energy level ei = eg + A/2 and sweep the other 
level, £2 = eo — A/2. This setup amounts to changing the 
field h externally, and is thus well suited for probing the 
magnetic response of our effective impurity. 

An example for such a protocol is presented in Fig. [Sj 
where ei is held fixed at ei — ~U/2. As £2 is swept 
through ci, a population inversion takes place, leading 
to a narrow dip in the conductance. The width of the 
conductance dip is exponentially small due to Kondo cor- 
relations. Indeed, one can estimate the dip width, Aedip, 
from the condition |ei — £2! — T/fle^^ei, which yields 



Aedip 



ttU_ 

' ST 



(58) 



B. Anisotropic couplings, Fi 7^ 

As demonstrated at length in Sec. IV Al the occur- 
rence of population inversion and a transmission zero for 
F| = F| requires an external modulation of the effec- 
tive magnetic field. Any practical device will inevitably 
involve, though, some tunnelling anisotropy, V| ^ V|. 
The latter provides a different route for changing the 
effective magnetic field, through the anisotropy-induced 
terms in Eq. 130(1 . Implementing the same protocol as 
in Sec. IV A 21 (that is, uniformly sweeping the average 
level position eg while keeping the difference A constant) 
would now generically result both in population inversion 
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FIG. 6: The occupation numbers (n;) and conductance G 
[in units of e^/{2nh)] as a function of eo + U/2 [in units of 
Ttot = (Fj -I- Fj)], calculated from Eqs. ilHt and lIH^ based 
on the mapping onto the Kondo model. The mod el p aram- 
eters are identical to those used in Fig. 2 of RefFLll lower 
left panel: h = ip = 0, (7/Ftot = 6, Ft/Ftot = 0.62415 and 
T = 0. The explicit tunnelling matrix elements are detailed in 
Eq. 15911 . corresponding to the rotation angles Oi — 2.1698 and 
9d ~ —0.63434 (measured in radians). The angle 6^ equals 
zero. The inset shows functional renormalization-group (fRG) 
data as defined in corrected for the renormalization 

of the two-particle vertex^S*— The small symbols in the in- 
set show the conductance as calculated from the fRG occu- 
pation numbers using our Eq. 1551 . The horizontal dotted 
lines in each plot mark the maximal conductance predicted 



by Eq. {e^ /2Tvh) ain^ 



and a transmission zero due to the rapid change in direc- 
tion of the total field htot- As emphasized in Sec. lIVB"2l 
the two phenomena will generally occur at different gate 
voltages when V-^ ^l- 



1. Degenerate levels, A = 6 = 

We begin our discussion with the case where A = 6 = 
0, which was extensively studied in RefUll It corre- 
sponds to a particular limit of the parallel-field config- 
uration where h = Q. In the parallel- field configuration, 
the conductance G and occupancies {rii) take the exact 
forms specified in Eqs. H44|) and (|54|l . respectively. These 
expressions reduce in the Kondo regime to Eqs. (|40(l and 
(|53|l . with dh either equal to zero or tt, depending on the 
sign of hl^^. 

Figure |H| shows the occupation numbers and the con- 
ductance obtained from Eqs. H4U|I and H53|) . for A = 6 = 
and the particular tunnelling matrix used in Fig. 2 of 
Refllj 



^(127 V(U6 



V(I33 -Va24 



(59) 



Here Aq equals \/Ftot/('?''p), with Ftot = T-f-t-Fj^ being the 
combined hybridization width. The Coulomb repulsion U 
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is set equal to GFtot , matching the value used in the lower 
left panel of Fig. 2 in Refflll For comparison, the cor- 
responding functional renormalization-group (fRG) data 
of Refll' is shown in the inset, after correcting for the 
renormalization of the two-particle vertexi^^'^^ The accu- 
racy of the fRG has been establishedii*i^ up to moderate 
values of C//Ftot ~ 10 through a comparison with Wil- 
son's numerical renormalization-group method: '^^ Includ- 
ing the renormalization of the two-particle vertex further 
improves the fRG data as compared to that of ReflTll. as 
reflected, e.g., in the improved position of the outer pair 
of conductance resonances. 

The agreement between our analytical approach and 
the fRG is evidently very good in the local-moment 
regime, despite the rather moderate value of U /Ftot 
used. Noticeable deviations develop in {rii) only as the 
mixed-valent regime is approached (for eq ^ ~Ftot or 
e-|-t/ < Ftot), where our approximations naturally break 
down. In particular, our approach accurately describes 
the phase lapse at eo = the inversion of popula- 

tion at the same gate voltage, the location and height of 
the correlation-induced resonances, and even the location 
and height of the outer pair of conductance resonances. 
Most importantly, our approach provides a coherent an- 
alytical picture for the physics underlying these various 
features, as elaborated below. 

Before proceeding to elucidate the underlying physics, 
we briefly quote the relevant parameters that appear 
in the conversion to the generalized Anderson model of 
Eq. Q. Using the prescriptions detailed in Appendix fXl 
the hybridization widths F^ — TrpV^ come out to be 

Ft/Ftot = 0.62415 , F jFtot = 0.36585 , (60) 
while the angles of rotation equal 



9i = 2.1698 , 



-0.63434. 



(61) 



Here 9i and 9^ are quoted in radians. Using the exact 
conductance formula of Eq. H44|) . G is predicted to be 
bounded by the maximal conductance 

2 2 

Gmax = 7^ sin^ 9i = 0.68210;^ , (62) 

obtained whenever the magnetization M — (n^ — ni)/2 
is equal to ±1/4. The heights of the fRG resonances are 
in excellent agreement with Eq. (|62|l . Indeed, as demon- 
strated in the inset to Fig. El the fRG occupancies and 
conductance comply to within extreme precision with the 
exact relation of Eq. (|55|l . As for the functional form of 
the Kondo temperature Tk, its exponential dependence 
on eo is very accurately described by Eq. In the 

absence of a precise expression for the pre-exponential 
factor when F| ^ Fj^ , we employ the expression 



Tk = iVUTtot/TT) 



exp 



7reo(t/ -I- ep) 
2C/(Ft - Fx) "'Fx 



in;^ 



(63) 



The occupancies and conductance of Fig.|Hlcan be fully 
understood from our general discussion in Sec. lIIII Both 
quantities follow from the magnetization M, which van- 
ishes at eo — —U/2 due to particle- hole symmetry. As 
a consequence, the two levels are equally populated at 
eo — —U/2 and the conductance vanishes [see Eqs. H44|l 
and H54(l ]. Thus, there is a simultaneous phase lapse and 
an inversion of population at eo = —U/2, which is a fea- 
ture generic to A = 6 = and arbitrary A. As soon as the 
gate voltage is removed from —U/2, i.e., eo = —U/2 + Se 
with (5e 7^ 0, a finite magnetization develops due to the 
appearance of a finite effective magnetic field ft-tot — hl^^z 
with 



F| - F I 1 + 25e/U 

■ In —'— 

1 - 25e/U 



(64) 



which properly reduces to Eq. ljB4|l (up to the small F^ 
correction in the exponent) when Ft = F| = F. 



[see Eq. (I33|) ]. Note that the sign of ft-tot coincides with 
that of (5e, hence M is positive (negative) for eo > —U/2 
(eo < —U/2). Since cos^d > for the model parameters 
used in Fig.|Bl it follows from Eq. (|54() that (rii) > (77,2) 
((ni) < (^2)) for eo > -U/2 (eo < -U/2), as is indeed 
found in Fig. |^ Once again, this result is generic to 
A = = 0, except for the sign of cos^^ which depends 
on details of the tunnelling matrix A. 

In contrast with the individual occupancies, the con- 
ductance G depends solely on the magnitude of M , and 
is therefore a symmetric function of 5e. Similar to the 
rich structure found for Ff = F| and A > in Fig. 01 the 
intricate conductance curve in Fig. |^ is the result of the 
interplay between /itot ^^.d Tk, and the nonmonotonic 
dependence of G on |M|. The basic physical picture is 
identical to that in Fig. ^ except for the fact that the 
effective magnetic field /itot now itself a function of the 
gate voltage eo. 

As a rule, the magnetization |M| first increases with 
\5e\ due to the rapid increase in ft-tot- It reaches its max- 
imal value Mmax at some intermediate \5e\ before de- 
creasing again as \5e\ is further increased. Inevitably |M| 
becomes small again once \5e\ exceeds U/2. The shape 
of the associated conductance curve depends crucially on 
the magnitude of Afmax, which monotonically increases 
as a function of U . When Mmax < 1/4, the conductance 
features two symmetric maxima, one on each side of the 
particle- hole symmetric point. Each of these peaks is 
analogous to the one found in Fig. 21 for A < he. Their 
height steadily grows with increasing U until the unitar- 
ity condition Mmax — 1/4 is met. This latter condition 
defines the critical repulsion Uc found in Ref^J- For 
U > Uc, the maximal magnetization M^ax exceeds one 
quarter. Hence the unitarity condition M = ±1/4 is 
met at two pairs of gate voltages, one pair of gate volt- 
ages on cither side of the particle-hole symmetric point 
eo — —U/2. Each of the single resonances for U < Uc 
therefore split in two, with the inner pair of peaks evolv- 
ing into the correlation-induced resonances of ReflU 
The point of maximal magnetization now shows up as 
a local minimum of the conductance, similar to the point 
eo = -U/2 in Fig. |3| when A > he. 
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For large U ^ Ftot, the magnetization \M\ grows 
rapidly as one departs from eg — —U /2, due to the expo- 
nential smallness of the Kondo temperature 7V|co=_c//2- 
The dot remains polarized throughout the local-moment 
regime, loosing its polarization only as charge fluctua- 
tions become strong. In this limit the inner pair of reso- 
nances lie exponentially close to eo — —U/2 (see below), 
while the outer pair of resonances approach |^e| ~ U/2 
(the regime of the conventional Coulomb blockade). 

The description of this regime can be made quantita- 
tive by estimating the position ±(5eciR of the correlation- 
induced resonances. Since M MK{h^^^^/Tj^) deep in 
the local-moment regime, and since (5eciR ^ Ttot for 
Ttot *C U, the correlation-induced resonances are peaked 
at the two gate voltages where h^^^ w ±2.4T^|(:g=_[//2. 
Expanding Eq. (|64|l to linear order in Sccm/U ^ 1 and 
using Ec^. (|63|) one flnds 



0.6- 



nU 



Ft -r 



-TK\eo = -U/2 



0.6— — exp 

i t — 1 



(65) 



functions at the Fermi energy: 



GAe = 0) = 



(67) 



Here, in writing the Green functions of Eq. H67|) . we 
have made use of the fact that the imaginary parts of 
the retarded dot self-energies, — F^, are unaffected by 
the Coulomb repulsion U at zero temperature at the 
Fermi energy. The energies e^r have the exact repre- 
sentations! ia = Fo- cot Scr in terms of the associated 
phase shifts 5„ — ■n{n„). Since M = ±1/4 at the peaks 
of the correlation- induced resonances, this implies that 
5a = 7r/2 ± (T7r/4, where we have set ritot — 1-^° Thus, 
the renormalized dot levels take the form = ^ctFo-, 
resulting in 



(68) 



The relation specified in Eq. H68|l was found in ReflU 
for the special case where Uli = a/si and aL2 = —o,r2^ 
Here it is seen to be a generic feature of the correlation- 
induced resonances for A = 6 = and arbitrary A. 



Here the pre-exponential factor in the final expression for 
(^eciR is of the same accuracy as that in Eq. (|63|l . 

We note in passing that the shape of the correlation- 
induced resonances and the intervening dip can be conve- 
niently parameterized in terms of the peak position (5eciR 
and the peak conductance G,nax- Expanding Eq. (|64|l to 
linear order in (5e/[/ <C 1 and using Eq. 1441) one obtains 



2. Nondegenerate levels: arbitrary A and b 



G{5e) = Gmax sin^ 



2tiM 



K 



2ASe 
SeciR 



(66) 



where MKih/Tx) is the universal magnetization curve of 
the Kondo model [given explicitly by This param- 

eterization in terms of two easily extractable parameters 
may prove useful for analyzing future experiments. 

It is instructive to compare Eq. for (Seem with the 
fRG results of ReflTH which tend to overestimate JeciR. 
For the special case where uli = uri and aL2 = —aR2, 
an analytic expression was derived for JeciR based on 
the fRG»ii The resulting expression, detailed in Eq. (4) of 
ReflU shows an exponential dependence nearly identical 
to that of Eq. (|65|l . but with an exponent that is smaller 
in magnitude by a factor of tt^/S « 1.23.^^ The same 
numerical factor appears to distinguish the fRG and the 
numerical renormalization-group data depicted in Fig. 3 
of Ref 11, supporting the accuracy of our Eq. (|65|l . It 
should be emphasized, however, that Fig. 3 of ReffTH 
pertains to the tunnelling matrix of Eq. (|59(l rather than 
the special case referred to above. 

We conclude the discussion of the case where A = b = 
with accurate results on the renormalized dot levels 
when the dot is tuned to the peaks of the correlation- 
induced resonances. The renormalized dot levels, e-\ and 
ii , can be defined through the T = retarded dot Green 



Once -f 6^ 7^ 0, the conductance and the partial 
occupancies can have a rather elaborate dependence on 
the gate voltage eo- As implied by the general discus- 
sion in Sec, mil the underlying physics remains driven by 
the competing effects of the polarizing field ft-tot and the 
Kondo temperature T^- However, the detailed depen- 
dencies on eo can be quite involving and not as revealing. 
For this reason we shall not seek a complete characteri- 
zation of the conductance G and the partial occupancies 
{rii) for arbitrary couplings. Rather, we shall focus on the 
case where no Aharonov-Bohm fluxes are present and ask 
two basic questions: (i) under what circumstances is the 
phenomenon of a phase lapse generic? (ii) under what 
circumstances is a population inversion generic? 

When ifL = ^R = 0, the conductance and the partial 
occupancies are given by Eqs. H40|l and H53|) . respectively. 
Focusing on G and on (ni — 77.2), these quantities share a 
common form, with factorized contributions of the mag- 
netization Mji and the rotation angles. The factors con- 
taining MK{htot/TK) never vanish when hsm.6 ^ 0, 
since /itot always remains positive. This distinguishes the 
generic case from the parallel-field configuration consid- 
ered above, where phase lapses and population inversions 
are synonymous with M — 0. Instead, the conditions for 
phase lapses and population inversions to occur become 
distinct once hsvuO 7^ 0, originating from the indepen- 
dent factors where the rotation angles appear. For a 
phase lapse to develop, the combined angle 9i+SR0h must 
equal an integer multiple of it. By contrast, the inver- 
sion of population requires that Od + sgOh — 7i'/2 modTr. 
Here the dependence on the gate voltage eo enters solely 
through the angle Oh, which specifies the orientation of 
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the effective magnetic field /itot [see Eq. H33|) ]. Since the 
rotation angles 9i and 9d are generally unrelated, this im- 
plies that the two phenomena will typically occur, if at 
all, at different gate voltages. 

For phase lapses and population inversions to be ubiq- 
uitous, the angle Oh must change considerably as eo 
is swept across the Coulomb-blockade valley. In other 
words, the effective magnetic field /itot must nearly flip its 
orientation in going from eg ~ to eq ~ ~U . Since the x 
component of the field is held fixed at h^^^ = hsvaO > 0, 
this means that its z component must vary from h^^^ 3> 
/iJq^ to — /ifot » /ifot as a function of eg- When this 
requirement is met, then both a phase lapse and an in- 
version of population are essentially guaranteed to occur. 
Since h^^^ crudely changes by 

A^^„t~-(rT-rjin(c//rtot) (69) 

TT 

as eo is swept across the Coulomb-blockade valley, this 
leaves us with the criterion 

(Ft - Fx) ln(C//Ftot) » . (70) 

Conversely, if IS? + P » (F| - Fx) ln([//Ftot), then 
neither a phase lapse nor an inversion of population will 
occur unless parameters are fine tuned. Thus, the larger 
U is, the more ubiquitous phase lapses become. ^'^^ 

Although the logarithm ln(J7/Ftot) can be made quite 
large, in reality we expect it to be a moderate factor of 
order one. Similarly, the difference in widths Ff — F^ 
is generally expected to be of comparable magnitude to 
F^. Under these circumstances, the criterion specified 
in Eq. lO reduces to F^ » y/A^TW. Namely, phase 
lapses and population inversions are generic as long as 
the (maximal) tunnelling rate exceeds the level spacing. 
This conclusion is in line with that of a recent numerical 
study of multi-level quantum dots4i 

Finally, we address the effect of nonzero h = y/A^+lP 
on the correlation-induced resonances. When h ^ 
F| ln(?7/Ftot), the effective magnetic field ft-tot ~ is 
large throughout the local-moment regime, always ex- 
ceeding F| and F|. Consequently, the dot is nearly fully 
polarized for all — C/ < eo < 0, and the correlation- 
induced resonances are washed out. Again, for practical 
values of U/Ttot this regime can equally be characterized 
by /i»F|ii 

The picture for F^ ln(t//Ftot) ^ h is far more elabo- 
rate. When TK\eQ=-u/2 3> h, the magnetic field is uni- 
formly small, and no significant modifications show up 
as compared with the case where h — 0. This leaves us 
with the regime TK\eo=-u/2 ^ ^ F|, where various 
behaviors can occur. Rather than presenting an exhaus- 
tive discussion of this limit, we settle with identifying 
certain generic features that apply when both compo- 
nents \hcos6\ and hsinO exceed TK\to=-u/2- To begin 
with, whatever remnants of the correlation-induced res- 
onances that are left, these are shifted away from the 
middle of the Coulomb-blockade valley in the direction 



where \h^^^\ acquires its minimal value. Consequently, 
htot and Tk no longer obtain their minimal values at the 
same gate voltage eo. This has the effect of generating 
highly asymmetric structures in place of the two sym- 
metric resonances that are found for ft, = 0. The heights 
of these features are governed by the "geometric" factors 
sin^(0; + sji9h) at the corresponding gate voltages. Their 
widths are controlled by the underlying Kondo tempera- 
tures, which can differ substantially in magnitude. Since 
the entire structure is shifted away from the middle of the 
Coulomb-blockade valley where Tk is minimal, all fea- 
tures are substantially broadened as compared with the 
correlation-induced resonances for h = 0. Indeed, similar 
tendencies are seen in Fig. 5 of RefO, even though the 
model parameters used in this figure lie on the borderline 
between the mixed- valent and the local-moment regimes. 



VI. CONCLUDING REMARKS 

We have presented a comprehensive investigation of 
the general two-level model for quantum-dot devices. A 
proper choice of the quantum-mechanical representation 
of the dot and the lead degrees of freedom reveals an ex- 
act mapping onto a generalized Anderson model. In the 
local-moment regime, the latter Hamiltonian is reduced 
to an anisotropic Kondo model with a tilted effective 
magnetic field. As the anisotropic Kondo model flows 
to the isotropic strong-coupling fixed point, this enables 
a unified description of all coupling regimes of the origi- 
nal model in terms of the universal magnetization curve 
of the conventional isotropic Kondo model, for which ex- 
act results are available. Various phenomena, such as 
phase lapses in the transmission phase^^iS charge oscilla- 
tions,^*— and correlation-induced resonances^^'^-^ in the 
conductance, can thus be accurately and coherently de- 
scribed within a single physical framework. 

The enormous reduction in the number of parameters 
in the system was made possible by the key observation 
that a general, possibly non-Hermitian tunnelling matrix 
A can always be diagonalized with the help of two si- 
multaneous unitary transformations, one pertaining the 
dot degrees of freedom, and the other applied to the lead 
electrons. This transformation, known as the singular- 
value decomposition, should have applications in other 
physical problems involving tunnelling or transfer matri- 
ces without any special underlying symmetries. 

As the two-level model for transport is quite general, 
it can potentially be realized in many different ways. As 
already noted in the main text, the model can be used to 
describe either a single two-level quantum dot or a double 
quantum dot where each dot harbors only a single level. 
Such realizations require that the spin degeneracy of the 
electrons will be lifted by an external magnetic field. Al- 
ternative realizations may directly involve the electron 
spin. For example, consider a single spinful level coupled 
to two ferromagnetic leads with non-collinear magneti- 
zations. Written in a basis with a particular ad hoc local 
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spin quantization axis, the Hamiltonian of such a system 
takes the general form of Eq. after properly combin- 
ing the electronic degrees of freedom in both leads. As is 
evident from our discussion, the local spin will therefore 
experience an effective magnetic field that is not aligned 
with either of the two magnetizations of the leads. This 
should be contrasted with the simpler configurations of 
parallel and antip arallel magnetizations, as considered, 
e.g., in R,efs ll7l2,^ andl27i 

Another appealing system for the experimental obser- 
vation of the subtle correlation effects discussed in the 
present paper is a carbon nanotube-based quantum dot. 
In such a device both charging energy and single-particle 
level spacing can be sufficiently large^ to provide a set of 
well-separated discrete electron states. Applying exter- 
nal magnetic field either perpendicular^^ or/and paral- 
lel^ to the nanotube gives great flexibility in tuning the 
energy level structure, and thus turns the system into a 
valuable testground for probing the Kondo physics ad- 
dressed in this study. 

Throughout this paper we confined ourselves to spin- 
less electrons, assuming that spin degeneracy has been 
lifted by an external magnetic field. Our mapping can 
equally be applied to spinful electrons by implementing 
an identical singular-value decomposition to each of the 
two spin orientations separately (assuming the tunnelling 
term is diagonal in and independent of the spin orienta- 
tion). Indeed, there has been considerable interest lately 
in spinful variants of the Hamiltonian of Eq. J^l, whether 
in connection with lateral quantum dots,^'^''*® capacitively 
coupled quantum dots^S'^S'^ or carbon nanotube de- 
vices^ Among the various phenomena that have been 
discussed in these contexts, let us mention SU(4) variants 
of the Kondo eSect^'^^^ and singlet-triplet transitions 
with two-stage screening on the triplet side.'*^''*^ 

Some of the effects that have been predicted for the 
spinful case were indeed observed in lateral semiconduc- 
tor quantum dots^ii^ and in carbon nanotube quantum 
dots44 Still, there remains a distinct gap between the 
idealized models that have been employed, in which sim- 
plified symmetries are often imposed on the tunnelling 
term, and the actual experimental systems that obviously 
lack these symmetries. Our mapping should provide a 
much needed bridge between the idealized models and 
the actual experimental systems. Similar to the present 
study, one may expect a single unified description encom- 
passing all coupling regimes in terms of just a few basic 
low-energy scales. This may provide valuable guidance 
for analyzing future experiments on such devices. 
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APPENDIX A: MAPPING PARAMETERS 

In this Appendix we give the details of the mapping 
of the original Hamiltonian, Eq. onto the generalized 
Anderson Hamiltonian of Eq. ([TJl . 

The first step is the diagonalization of the matrix A, 
Eq. which describes the couphng between the dot 
and the leads in the original model. Since A is generally 
complex and of no particular symmetry, it cannot be di- 
agonalized by a single similarity transformation. Rather, 
two (generally different) unitary matrices, Rd and Ri, are 
required to achieve a diagonal form. 



Ft 
Vi 



RiAR\ 



(Al) 



This representation, known as the singular-value de- 
composition, is a standard routine in numerical pack- 
ages. Here we provide a fully analytical treatment of the 
2x2 case relevant to our discussion. To this end we 
parametrize the two rotation matrices in the form 



Rd = e'^^e'^^-/^'>^^U{ed,cj)d), 



where 



cos(6i/2) e-"^sin(6l/2) 
e^'^sin(6'/2) cos(6l/2) 



(A2) 
(A3) 



(A4) 



describes a rotation by angle 9 about the axis — sin((/)) x+ 
cos(0) y. 

The various parameters that enter Eqs. HA2|I and (|A3() 
have simple geometrical interpretations. The two sets 
of angles, {Od,4>d) and {0i,4>i), are the longitudinal and 
the azimuthal angles of the vectors pointing along the 
direction of the z axis which defines the corresponding 
spin variables in Eq. iQ, see Fig. [7| The three angles 
XA, Xd, and xi correspond to the choice of the phases 
of the single-particle operators dj. and c\^. The lat- 
ter angles are chosen such that the matrix elements of 
the transformed Hamiltonian, Eq. (jT)), will be real with 
hsvaO > 0. Note that Rd and Ri are determined up 
to a common overall phase. This degree of freedom has 
been exhausted in Eqs. ljA2|l and ljA3|l by requiring that 
det Ri = 1. 

In order to determine the rotation matrices Rd and Ri , 
one diagonalizes the hermitian matrices AA'' and A'' A, 
whose eigenvalues are evidently real and equal to |VJjp. 
This calculation determines the matrices U{9d,4>d) and 
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where 



FIG. 7: The original dot degrees of freedom, d\ and d^, define 
a pseudo-spin- i representation with the axes x' , y' , and 2'. 
The level indices 1 and 2 are identified in this representation 
with ±i spin projections on the z' axis. The energy splitting 
A and the hopping b combine to define the magnetic-field 
vector hm . The unitary matrix Rd takes the spin to a new 
coordinate system whose axes are labelled with x, y, and z. 
The new "spin" labels a =t and a =1 represent ±i spin 
projections on the new z axis, whose direction is defined by 
the longitudinal and the azimuthal angles 9d and <j)d- The new 
X axis lies in the plane of vectors z and hm' . A similar picture 
applies to the conduction-electron degrees of freedom, where 
the lead index a — L, R plays the same role as the original 
level index i = 1, 2. 



U{9i^(f)i), and yields the values of |Vo-|. Indeed, using 
Eqs. (HH), ijUll, and l|X3)l . one obtains 



= U{edAd)A^AU\ed,(j>d). (A5) 



Assuming > |Vx| (the case where \Vt^\ — \Vi\ is 
treated separately in Sec. I A 2p . these two equations give 



(A6) 



\VJ'=X±Y . 



2 arctan . 



d/l 



Y + Z, 



d/l 



(A7) 



arctan ( "^^"^^ °-RiQfl2 'f_ 1 ^ ^ ^ ^^g-j 

\aLiaL2 + aRiaB2 2 



oi — arctan | ; tan — ■ I + tt 77; , (A9j 



\aL2aR2 + aLio-Bi 2 



Y = 
Zi = 



X2- |detA|2, 



2) ) 



a=L,R 



2 

Li 



(AlO) 

(All) 
(A12) 

(A13) 



and 



{aLiaL2 + aRiaB2) cos — , (A14) 
(aiiflfli +aL2afl2)cos^ . (A15) 



2r]d = 1 - sgn 
2-qi = 1 - sgn 



The plus sign in Eq. (|A6p corresponds to V|, since the 
spin- up direction is defined as the one with the larger cou- 
pling, |V^||2 > \Vi\'^. The longitudinal angles < 6*^,6'; < 
TT are uniquely defined by Eq. (|A7I) . while the quadrants 
for the azimuthal angles — 7r/2 < (fid, (f'l ^ 37r/2 must be 
chosen according to Eqs. (|A14p and IjAlSp . The auxil- 
iary quantities in Eqs. (jA10p - (jA13p obey the inequalities 
X>Y i,ndY>\Zd/i\. 

The next step is to determine the angles XA, Xd, and xi 
which come to assure, among other things, that V| > 
are both real and non- negative. Let us begin with XA- 
When det A ^ 0, i.e., for V| > 0, the angle XA is uniquely 
determined by taking the determinants of both sides of 
Eq. IjAljl and equating their arguments. This yields 



XA = ^ argdeti. 



(A16) 



When det A — 0, the angle XA can take arbitrary values. 
This stems from the fact that Vi — 0, and therefore c^i 
can be attached an arbitrary phase without affecting the 
form of Eq. (|7)|. In this case we choose XA = 0. 

Next we rotate the Hamiltonian term £d, which is the 
first term of the isolated dot Hamiltonian, Eq. (j2Jl . Upon 
converting to the rotated dot operators and dj, the 
single-particle term £d transforms according to 



£d 



Rd ^d ^d 



(A17) 



Consider first the partial rotation U{9d, <j>d)£dU^ {Od, 4>d) 
[see Eq. (IA2|) ]. Writing £d [as defined in Eq. igj] in the 
form 



£d — £0 



TO • T 



with 



and 



TO 



h = v/A2 + 52 



(A18) 
(A19) 



-cos x-h-T-sm y--rz, (A20) 
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the partial rotation U{9d,(t)d)^dU^9d,4'd) gi 



ives 



eo - - m • T , 



(A21) 



where m is the unit vector obtained by rotating m' by an 
angle —9d about the axis — sin((/)d) x+cos{4>d) V- Defining 
the angle 9 G [0, tt] which appears in Eq. lO according 
to cos 9 ~ rriz , it follows from simple geometry that 

A 

cos 9 = — — cos 9d 
li 

- ^ sm9d cos[(j)d + ifL - </'fl)/2] • (A22) 
n 

The full transformation R^EdR\ corresponds to yet an- 
other rotation of ra by an angle — Xd about the z axis. 
The angle Xd is chosen such that the projection of rh 
onto the xy plane is brought to coincide with the x di- 
rection. This fixes Xd uniquely, unless hsa\9 happens to 
be zero (whether because /i = or because 9 is an integer 
multiple of tt). When hsm9 = 0, the angle Xd can take 
arbitrary values. Physically this stems from the fact that 
spin-up and spin-down degrees of freedom can be gauged 
separately within Eq. (0 . We choose Xd = in this case. 
The explicit expression for Xd when h sin ^ is quite 
cumbersome and will not be specified. As for the remain- 
ing angle xu it is fixed by the requirement that Va will 
be real and non-negative. 

Note that the conditions for the two exactly solvable 
cases quoted in the main text, Eqs. ((Tn|l and l(TT)) . are 
readily derived from our expressions for the eigenvalues 
Va. The first case, Eq. l(TU)l. corresponds to — 0, 
which requires det A = e'^'^aLiaii2 — aL2a/si = 0. This 
immediately leads to Eq. (|10() . The second solvable case, 
Eq. 111() . corresponds to equal eigenvalues, which implies 
y = [Eqs. (|X6)) . I|A10|) and I|A1HI remain intact for 
\V-\ \ = By virtue of the inequalities Y > \Z^ii\, this 

necessitates that Zd and Z; are both zero, which gives 
rise to the first two conditions in Eq. (|ll|l . The remaining 
condition on the Aharonov-Bohm phase follows from 
substituting the first two conditions into the definition of 
Y and equating Y to zero. 



1. No Aharonov-Bohm fluxes 

Of particular interest is the case where no Aharonov- 
Bohm fluxes are present, lpl — fR — 0. In the absence 
of a real magnetic field that penetrates the structure, the 
parameters that appear in the Hamiltonian of Eq. Q are 
all real. This greatly simplifies the resulting expressions 
for the rotation matrices Rd and Ri, as well as for the 
model parameters that appear in Eq. Q. In this subsec- 
tion, we provide explicit expression for these quantities 
in the absence of Aharonov-Bohm fluxes, focusing on the 
case where > Vi- The case where V-^ = Vi is treated 
separately in Sec. I A 21 



As is evident from Eqs. (jA8|) and ljA9|l . each of the az- 
imuthal angles 0^ and (f>i is either equal to or tt when 
(fi ~ 0. (The corresponding y' and y axes are parallel 
in Fig. [7|) It is therefore advantageous to set both az- 
imuthal angles to zero at the expense of extending the 
range for the longitudinal angles 9d and 9i from [0, tt] 
to (— 7r,7r]. Within this convention, Eq. ljA7p is replaced 
with 



V' = 2sd/ia,rcta,n, 



Y 


- Zd/i 




+ Zd/i 



where 



Sd = sgn{aLiaL2 + aRiaR2) 
si = sgn {a Lia Ri + a L2aR2) 
Similarly, the unit vector rh' of Eq. (IA2 

A 



(A23) 



(A24) 
(A25) 

reduces to 



m = X z 

h h 



(A26) 



which results in 



■ cos 9d 



■ smt 



b . 

— sml 
h 



— cos 9d 
n 



and 



arccos 



A 



■ sm Ud 



■ cos 9d 



(A27) 



(A28) 



Since the rotated unit vector m has no y component, 
the angle Xd is either equal to or tt, depending on the 
sign of m^. Assuming det A ^ and using Eqs. I)A4|I and 
(|A16|I . one can write Eq. (|A2I) in the form 



Rd (sgn det A) ' e 



1/2 ,^(l-se)r,/4 J(ea/2)7 



(A29) 



Se = sgn rux = sgn (A sin 9d ~ b cos ( 



(A30) 

Note that the flrst exponent in Eq. (|A29p is equal to 1 
for sg = +1, and is equal to e^^'^^'^^'^^ for sg — —1. If 
det^ = we set sgn det A 1 in Eq. HA29I) . while for 
A sin 9 d = bcos9d we select sg ~ +1. 

Proceeding to the remaining angle xi > we note that Rd 
of Eq. ljA29(l is either purely real or purely imaginary, 
depending on whether 



SR — Sg sgn det A 



(A31) 



is positive or negative. Since both e*'^^'/^-''^'' and A are 
real matrices, then e*'^-^'/^-''^^ must also be either purely 
real or purely imaginary in tandem with Rd in order for 
the eigenvalues and Vi to be real. This consideration 
dictates that xi is an integer multiple of tt, with an even 
(odd) integer for positive (negative) sr. The end result 
for Ri is therefore 



Ri 



(A32) 



Here ry^ = ±1 is an overall phase which comes to assure 
that the eigenvalues Vij and Vi are non-negative. 
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2. Isotropic couplings, V| = Vi 

Our general construction of the rotation matrices Rd 
and Ri fails when = \Vi\ — V. Equations HA6|I . 
(|A10|) and (|A11|I remain in tact for \V^ \ = | V| | , however 
the angles d^/i and (pd/i are ill-defined in Eqs. ljA7p - ljA9p . 
This reflects the fact that the matrices A^A and AA'' are 
both equal to V'^ times the unit matrix, hence any rota- 
tion matrix U{6,(j)) can be used to "diagonalize" them. 
There are two alternatives for treating the isotropic case 
where |V|| = |V||. The first possibility is to add an in- 
finitesimal matrix ijB that lifts the degeneracy of \V^\ 
and \Vi\: A A + ijB. Using the general construc- 
tion outlined above and implementing the limit 77 — > 0, a 
proper pair of rotation matrices Rd and Ri are obtained. 
The other alternative is to directly construct the rota- 
tion matrices Rd and Ri pertaining to this case. Below 
we present this second alternative. 

A key observation for the isotropic case pertains to the 
"reduced" matrix 



f= (deti) ^^^i, 



which obeys 



T^T = TT^ = 1 



det T = 1 . 



(A33) 



(A34) 



As a member of the SU(2) group, T can be written in 
the form 



T = U{eTAT)e 



i(XTl2)i 



(A35) 



with Ot G [0, tt]. Explicitly, the angles 9t, (f>T and xt are 
given by 

Ot = 2arccos(|detir^/^ |aLi|) , (A36) 
XT = 2 arg[(det A)-'^^ gli] , (A37) 



and 



br = arg[(det A)-'^^ am] - vr - xt/2 . (A38) 



Exploiting the fact that det A — V'^, the matrix A takes 
then the form 



A^V R\Rd = V e'^^ U{0T, cj)T) e 



i(XT/2)T 



(A39) 



where the angle XA is defined in Eq. IjAlBjl . 

Equation ljA39p determines the matrix product RjRd- 
Any two rotation matrices that satisfy the right-most 
equality in Eq. ljA39p transform the tunnelling matrix 
A to V times the unit matrix, as is required. The rota- 
tion matrix Rd is subject to yet another constraint, which 
stems from the requirement that hsinO > in Eq. ((TJ. 
We note that this constraint as well does not uniquely 
determine the matrix Rd^ Perhaps the simplest choice 
for Rd is given by 



with 

1 TT 

Xd = -^{'fL - ^r) + -{1 ~ sgnb) , 

which corresponds to 

hcos8 ^ -A , hsmO |&| . 



(A41) 



(A42) 

Adopting the choice of Eq. ljA40|) , the rotation matrix Ri 
takes the form 

where 0t, 4>t, Xt and Xd are listed above. 

APPENDIX B: BETHE ANSATZ FORMULAE 

In this appendix we gather for convenience all relevant 
details of the exact Bethe ansatz solutions for the impu- 
rity magnetization in the isotropic Kondo and Anderson 
models in the presence of a finite magnetic field. Exten- 
sive reviews of these solutions (including the anisotropic 
Kondo model) are available in the literature. ^^•^''■^^ Here 
we only summarize the main results of relevance to our 
analysis, and briefly comment on the numerical proce- 
dure. We confine ourselves to zero temperature, al- 
though explicit equations do exist also at finite temper- 
ature. Throughout the Appendix we employ units in 
which 1139= 1- 

1. Isotropic Kondo model 

We begin the presentation with the case of a Kondo 
impurity, before turning to the more elaborate case of an 
Anderson impurity. As a function of the magnetic field /i, 
the magnetization of an isotropic spin-^ Kondo impurity 
is given by the explicit expression [see, e.g., Eq. (6.23) of 
ReflB 



M{h) 



40F 



+°|^ {iuj + QY'^/'^'' sech(f ) ( 



(Bl) 



Here T{z) is the complex gamma function. The Kondo 
temperature, Tr-, is defined via the inverse of the spin 
susceptibility. 



T^^ = 27r lim M{h)/h 



(B2) 



Evidently, Eq. ||BTJ is a universal function of the ratio 
h/Tx, which is denoted in the main text by Mxih/TK)- 
It has the asymptotic expansion 



M{h) 



h/{2TTTK) 



A\n(h/TH) 



lnln(/i/Tff ) 
81n^(h/T„) 



(B3) 



Rd = e'^^e 



(A40) where Th 
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2. Isotropic Anderson model 

In contrast to the Kondo model, there are no closed- 
form expressions for the total impurity occupancy n^ot = 
(n| +n^) and magnetization M = (n| — nj^)/2 in the 
isotropic Anderson model. The exact Bethe anstaz solu- 
tion of the model provides a set of coupled linear integral 
equations from which ritot and M can be computed. Be- 
low we summarize the equations involved and comment 
on the numerical procedure that is required for solving 
these equations. The expressions detailed below apply to 
arbitrary eo, U, h, and F at zero temperature. In terms 
of the Hamiltonian of Eq. (TJ , we restrict the discussion 
to F = Ff = Fj (isotropic Anderson model) and sin 9 — 
(parallel- field configuration). The case sinO ^ follows 
straightforwardly from a simple rotation of the dot and 
the conduction-electron operators about the y axis. 

a. Kondo temperature 

The most accurate analytical expression that is avail- 
able for the Kondo temperature of the isotropic Anderson 
model can be written as 

Tk = (\/2[7f/^) exp[7r (F^ + eq t/ + e') /(2 t^T)] , 

(B4) 

where F — irplVl"^. This exp ression for Tk exactly repro- 
duces Eq. (6.22) of Ref[23 for the symmetric Anderson 
model, eo = —U/2. It also coincides with Eq. (7.11) 
of Ref 22 for the Kondo temperature of the asymmetric 
model when U ^ T. Note that the F^ term in the expo- 
nent is usually omitted from Eq. ljB4|) on the basis of it 
being small. It does in general improve the estimate for 
Tk- 

In the local-moment regime, where Eq. (|B4|I is valid, 
the impurity magnetization of the isotropic Anderson 
model is dominated by the universal magnetization curve 
of Eq. (|Bll) up to fields of the order of - VW > Tk 
(see, e.g., lower left inset to Fig.|21). At yet larger fields, 
h 3> VTU, the magnetization of the Anderson model can 



no longer be described by that of the Kondo model, as 
charge fluctuations become exceedingly more important 
than spin flips. Rather, M is well described by perturba- 
tion theory in F. Importantly, the asymptotic expansion 
of Eq. ljB3|) properly matches (to leading order in T/U) 
the perturbative resul1»2i for M when h ~ \/TU. Thus, 
the two approaches combine to cover the entire range in 
h for the Anderson model. 



b. Bethe ansatz equations for the occupancy and 
magnetization 

The Bethe ansatz solution of the Anderson model fea- 
tures four key quantities, which are the distributions of 
the charge and spin rapidities, Pi/h(fc) and a'i/h(A), re- 
spectively, for the impurity (i) and the host (h) band. 
The total impurity occupancy and magnetization are ex- 
pressed as integrals over the distributions of the charge 
and spin rapidities for the impurity: 

M ^ I [ Pi{k)dk, (B5) 
rQ 

Ud = I- ai{X)dX. (B6) 

J —OO 

The upper limits of integration in Eqs. (|B5p and ljB6p 
are determined through implicit conditions on the corre- 
sponding distribution functions for the host band, 

A ^ j"" ~p^{k)dk, (B7) 

= r d^{\)d\. (B8) 

As for the distributions of the rapidities for the impu- 
rity and the host, these are determined by the same pair 
of linear integral equations, only with different inhomo- 
geneous parts: 



where^ 

S{x) 

R{x) 
9{k) 



dgjk) 
dk 



^.,.nJ^.,J , Mk) ["^ 



R[g{k)- g{k')]p{k')dk' + 



dk 



S[g{k)-X]aiX)d\^p<^'>Hk), 



a{X)- I R{X~ X')a{X')dX' + I S[X~g{k)]p{k)dk^a^°\X), 

I 



(B9) 
(BIO) 



(here ^' is the digamma function). The inhomogeneous 



1 



2 cosh(7ra;) 
1 











^ (l 


+^2) 















= — Re 

^ {k-ep- U/2f 
2UV 



(Bll) 
(B12) 
(B13) 
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parts in Eqs. (|B9|I and (jB10|) are given in turn by 

P^\k) = A(fc) + ^ £y^9ik) - g{k')] A(fc') dk' , 

(B14) 

(B15) 

/ + 00 
S[\ ~ g{k)] A{k) dk , (B16) 
-OO 

4''^(A) = 2^ J J[X-9{k)]dk, (B17) 
where A(fc) is the Lorentzian function 

A{k) = -— -. (B18) 

c. Details of the numerical procedure 

The main obstacle faced with in a numerical solution 
of the Bethe ansatz equations is the self-consistent de- 
termination of the upper integration bounds that appear 
in Eqs. ljB5|l - ljB10p . These are computed iteratively ac- 
cording to the scheme 

pI"~'^ , 4""'^ => B^-^ , Q(") ^ pI"^ , . (B19) 

Starting with p'^ '^^ and ct^" as input for the nth it- 
eration, and Q*^"-' are extracted from Eqs. (|B7|1 and 



(|B8|) . Using the updated values for B and Q, p]^ (k) and 

4""* (A) ^^'^ then obtained from the solution of Eqs. ljB9|) 
and (|B10I) . This cycle is repeated until convergence is 
reached. The first iteration in this procedure is usually 
initialized with p'^^k) and ^''''(A) as input. Standard 
techniques are then used to ensure rapid convergence of 
the iterative solution. Typically 15 to 30 iterations are 
required to achieve a relative accuracy of 10~^ for the 
vector {B,Q). 

The core of this cycle is the solution of Eqs. I|B9I) and 
(|B10|) . These are solved (for given values of B and Q) 
by discretizing the integration interval with adaptively 
chosen 500 1000 mesh points. Once a self-consistent 
solution is reached for B, Q, ph{k) and (Th(A), the corre- 
sponding distributions of rapidities for the impurity are 
obtained from a single solution of Eqs. I|B9|I and l)B10|l . 
The impurity occupancy and magnetization are calcu- 
lated in turn from Eqs. HB5|) and HB6I) . 

To test the accuracy of our numerical results, we have 
extensively checked them against the analytical solution 
for the zero-field occupancy nd{h — 0) and the zero- field 
susceptibility dM/ dh\fi=o. In suitable parameter regimes, 
we have also compared our results to perturbation the- 
ory in both U and F. In all cases tested the relative 
errors in rid and M were less than 0.05% and 0.5%, re- 
spectively. This accuracy can be systematically improved 
by increasing the number of discretization points used in 
solving Eqs. I|B9|) and IjBlOjl for the distributions. Our 
results were also in full agreement with those reported by 
Okiji and Kawasaki, except for M{h) where up to 10% 
differences were found. Considering the extensive set of 
checks that were applied to our results, it appears that 
the discrepancy is due to lower numerical accuracy in the 
solution of Ref. i23i 
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